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The muhiphcity, average transverse momentum, and charged particle 
transverse momentum distributions have recently been measured in LHC 
experiments. The multiplicity and average transverse momentum grow 
with beam energy. Such growth is expected in the theory of the Color Glass 
Condensate, a theory that incorporates the physics of saturation into the 
evolution of the gluon distribution. We show that the energy dependence 
of the pp data and the LHC data for pp scattering at y/s > 200 GeV 
may be simply described using a minimal amount of model input. Such 
a description uses parameters consistent with the Color Glass Condensate 
descriptions of HERA and RHIC experimental data. 



1. Introduction 

The first LHC data has been released on total charged particle multiplic- 
ity as a function of energy, average transverse momenta of charged particles 
as a function of energy p^-fS] and transverse momentum of charged particles 
as a function of charged particle multiplicity [6] . The generic features of the 
experimental data are that the charged particle pseudo-rapidity densities 
rise as a power of energy, dN/drj ~ E^*^-^^, and that the average transverse 
momentum rises with both energy and charged particle multiplicity. This 
behavior has a natural explanation within the theory of saturation and the 
Color Glass Condensate |7j -tl2j- Within this theory, the total multiplicity 
of produced particles is computable [13]- [H]- Assuming local parton hadron 
duality |T2] , the initially produced gluon multiplicity is proportional to the 
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charged particle multiplicity. The pseudo-rapidity multiplicity density can 
be expressed in terms of the saturation momentum as |16j-|21j: 

1 dNch _ const. 2 (-.^ 

The strong coupling constant is evaluated at the saturation momentum 
scale. 

The saturation momentum scale is proportional to the transverse gluon 
density. Evolution equations that build in the effects of high gluon density 
have been derived in Refs. [22j-[33j. The dependence of the saturation 
momentum on fractional gluon momentum may be determined from such 
considerations [3l]. It has been shown that the generic features of such 
a description of saturation can describe the HERA data on inclusive and 
diffractive deep inelastic scattering [35j-t38j, and correctly predicts observed 
scaling properties of experimental data [39] -[lO]. The results of this analysis 
relevant for our purposes is that the saturation momentum at x values 
appropriate for the LHC scales with x as Q2^t ~ l/x^ where A ~ 0.2 - 0.3 

There have been detailed saturation based predictions and descriptions 
of the recent results from the LHC [21], [1I]-[13]- The goal of this paper is 
not to improve upon the descriptions provided in these works. Our goal is 
to show that the simplest generic features of saturation based descriptions 
are adequate to quantitatively describe data on the dependence particle 
multiplicities on energy and the dependence of such average transverse mo- 
mentum upon multiplicity and beam energy. We will also argue that there 
is an approximate geometric scaling of transverse momenta distributions 
measured at LHC energies. 



2. Color Glass Condensate Description of the LHC data 

To reduce the Color Glass Condensate description to its simplest possible 
form, we will assume that the density of produced charged particles per unit 
pseudo-rapidity scales as 

^ = .Ql, = AE' (2) 

Here the energy is measured in units of TeV. 

The parameters A and A can be determined by a fit to the LHC data. In 
Fig. [H an excellent fit is shown to the data that also includes lower energy 
data for proton-anti-proton scattering. (This agreement with the pp data 
is a little surprising since there should be some small difference between pp 
and pp scattering at the energies of interest.) 
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Fig. 1. Charged particle rapidity density as a function of energy compared to a 
power law 3.7043 E^-^^ Data from the LHC[T]-[5], UA(1)[11], UA(5)g5] and CDF 

[ISl-gT]. 



On dimensional grounds, at very high energy, we expect the average 
transverse momentum will be proportional to the saturation momentum. 
Since this term does not entirely dominate the contribution to the transverse 
momentum at accessible energies, we add a constant, so that the result has a 
reasonable low limit at lower energy. In this case, there is a small difference 
seen in the data between pp and pp scattering. The power law growth should 
be, however, universal. We fit the average pT with the following form 



(pt) = B + CE^/'^. 



(3) 



The results of such a fit are shown in Fig. [21 We see that indeed the power 
of the energy and its coefficient (within errors which we do not quote here) 
are identical both for pp and pp data. Similar form of the average (px) has 
been recently postulated in Ref. [35] with higher power of E equal to 0.414. 
However, the recent 7 TeV CMS point is far below their curve. 

If the saturation momentum is the only scale that controls p^ distribu- 
tions, on dimensional grounds, these distributions should have a geometrical 
scaling 

1 dNch ^ f PT \ 



(5) 



- P 

a drid?pT: \ Q^at {pt / y/s) , 
This means that there is a universal function of the scaling variable 
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Fig. 2. Average transverse momentum as a function of energy compared to a power 
law B + C£;" "5 Data from the LHC[T]-[S], UA(1)[11], UA(5)[1S] and CDF gS]- 



that describes the data at different energies (where and ^/s are in GeV) . 
This is of course limited in the range of so that one is not probing quark 
and gluon distributions outside the saturation region. By a rescaling of 
variables, one can check if the data from CMS falls on a universal scaling 
curve. Since the data points of the CMS pT distributions are not yet publicly 
accessible we shall use throughout this paper an analytical parametrization 
in terms of Tsallis fit [35] as given in Refs. [U [S]: 

dN,^ = (l + ^\ ~" (6) 

drjcPp'Y E drj \ nT J 



where -Et = a/"^^ +p\ — "t-tt and 
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(7) 



Up to about 4 — 6 GeV, the limit of the available data, such a scaling 
relation is satisfied, as shown in Fig. O 

The scaling plot in Fig. [3] has been obtained by using power A in Eq.([2]) 
fixed from the DIS data at HERA. It is interesting to see whether this is also 
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Fig. 3. The CMS data for transverse momentum distributions on the left as func- 
tions of p\ (note that for A = scahng variable r = p\/{l GeV^ ) ) . On the right , the 
same px distribution rescaled in terms of the scaling variable r = Px/Qsat (-Pt / \/s) ■ 



the optimal power for hadron-hadron scattering. To this end we compute 
the mean deviation of the scaled distributions for different energies 



1 diVch 
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and normalizing them to the sum 



SE1-E2 




dijdr 



dN^h 



we define quantities 



A 



E1-E2 



crEi-E2 



dr 



E2, 



(9) 



(10) 



SE1-E2 

that are plotted in Fig. HI We see that the minima obtained with the Tsallis 
fit d^D are rather shallow and include the optimal value of A obtained from 
DIS, although the preferred value would be slightly bigger. We checked, 
however, that the higher value of A is incompatible with the energy depen- 
dence of charged multiplicity shown in Fig. [TJ We find this agreement (note 
that we use Tsallis parametrization instead of real data) as a strong support 
of the applicability of geometric scaling to hadron-hadron scattering. 

In the ATLAS experiment, the average of events with various mul- 
tiplicities was computed. There was a transverse momentum cutoff of 
Pt > 500 MeV. We expect as in Eq. ([3]), that the average transverse 
momentum will be 

{PT) = C + Dy^ (11) 
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Fig. 4. Normalized square deviations between scaled pt distributions for different 
CMS energies as function of saturation parameter A. Optimal DIS A = 0.23 is 
marked by a thin vertical line. 



Since the CMS and ATLAS cuts are different we simply show in Fig. [5] 
that such functional form provides a good description of the experimental 
data. We have presented two fits: one to the whole region of available 
multiplicities and the second one for N^h > 15. The latter choice is dictated 
by the slight change of the curvature of the data around N^h ~ 10. 




Fig. 5. ATLAS data compared to the square root fit of Ea. ((TT|) . Short-dashed 
red curve corresponds to the fit over the whole range of multiplicities, whereas 
long-dashed blue fit is restricted to A^ch > 15. 
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Fig. 6. CMS particle spectra in terms of Tsallis paranietrization of Refs. [H [5] for 
different energies. Solid line and long and short dashed lines below correspond to 
"v/i = 7, 2.36 and 0.9 TeV respectively. Dashed and short dashed lines above the 
solid hne correspond to our predictions for y/s = 10 and 14 TeV respectively. 

3. Predictions for Higher Energy 

Using the scaling analysis in this paper we can make predictions for the 
multiplicity per unit rapidity, average transverse momentum and transverse 
momentum distributions at higher LHC energies. Our predictions for the 
multiplicity per unit rapidity are shown in fact in Fig.[TJ In order to estimate 
roughly the error of that fit we simply propagated the experimental error of 
the 7 TeV point with the help of Eq.© obtaining dNch/drj\^^Q = 6.29±0.25 
and 6.80 it 0.27 for ^/s = 10 and 14 TeV respectively. In a similar way 
we have estimated average (px) = 0.562 zt 0.017 and 0.579 ± 0.017. By 
minimizing A10-7 and Ai4_7 with respect to parameters T and n of the 
Tsallis formula ([6]) we have obtained pT distributions at y/s = 10 and 14 TeV 
that are shown in Fig. The corresponding Tsallis parameters read: Tiq = 
0.153, mo = 6.6 for ^ = 10 TeV and Tu = 0.162, riu = 6.7 for 14 TeV. 

The scaling behavior we see in pp collisions can be used to estimate 
initial state effects for the heavy ion collisions. Such effects might be very 
important for measuring jet quenching effects, once the A dependence of 
the saturation momentum is established at LHC energy. Note that the 
asymptotic behavior of the Tsallis fit for high energies is controlled by the 
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parameter n^T/pT- If scales as A^^^, then the asymptotic hmit is obtained 
only at very high transverse momentum values, suggesting that saturation 
effects can influence transverse momentum distributions out to very high 
large values. This may influence experimental studies of jet quenching as 
an attempt to extract properties of the Quark Gluon Plasma. 
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